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Abstract 



We review the uncertainties in the spin-independent and -dependent elastic scat- 
tering cross sections of supersymmetric dark matter particles on protons and neu- 
trons. We propagate the uncertainties in quark masses and hadronic matrix elements 
that are related to the vr-nucleon a term and the spin content of the nucleon. By far 
the largest single uncertainty is that in spin-independent scattering induced by our 
ignorance of the {N\qq\N) matrix elements linked to the vr-nucleon a term, which 
affects the ratio of cross sections on proton and neutron targets as well as their abso- 
lute values. This uncertainty is already impacting the interpretations of experimental 
searches for cold dark matter. We plead for an experimental campaign to determine 
better the n-nucleon a term. Uncertainties in the spin content of the proton affect 
significantly, but less strongly, the calculation of rates used in indirect searches. 
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I. INTRODUCTION 



The most convincing way to confirm the existence and nature of dark matter parti- 
cles would be to observe directly their scattering on nuclei in low-background underground 
experiments. The sensitivities of these experiments are currently improving rapidly, and 
beginning to cut into the parameter space of plausible supersymmetric scenarios In 
order to evaluate accurately the impacts of these experiments, it is important to understand 
and minimize the hadronic uncertainties in the elastic scattering matrix elements for any 
given supersymmetric model. The rates for elastic scattering also control the rates for the 
capture of dark matter particles by celestial bodies such the Sun or Earth. There are good 
prospects for increasing significantly the sensitivities of experiments looking indirectly for 
astrophysical dark matter via the products of their annihilations in such bodies, adding to 
the motivations for understanding and reducing their uncertainties. 

Beyond the interpretation of upper limits on dark matter scattering may lie the inter- 
pretation of any eventual detection of a signal and the task of identifying the nature of 
the dark matter particle. In principle, there are four observables that could contribute to 
such an analysis, namely the spin-independent and -dependent cross sections on protons 
and neutrons, respectively. Part of the strategy for identifying the nature of any detected 
dark matter signal would be the comparison of the measured rates for scattering on different 
targets, with the comparison of spin-independent and -dependent scattering rates playing a 
particularly important role as a diagnostic tool [2I, S] • As we see later, currently there are 
considerable uncertainties also in such comparisons, related principally to uncertainties in 
the hadronic matrix elements of higher- dimensional effective interactions. 

In this paper, we consider only hadronic uncertainties in the elastic scattering rates. 
There are also potentially important uncertainties related to the supersymmetric model 
itself, namely how accurately one can estimate the coefficient of a given higher-dimensional 
effective interaction in a given model, and also in the astrophysical density of dark matter 
particles. Reducing the model uncertainty would require, e.g., a complete calculation of 
radiative corrections to the effective scattering operator, which lies beyond the scope of this 
work. As for the local density of cold dark matter, it is usually taken to be 0.3 GeV/cm^, 
but lower values have occasionally been advocated. 

The hadronic uncertainties we consider are listed in Table [11 They include those in the 
quark masses, expressed as md c,b,t and the ratios m^/m^ and ms/rrid, those in the matrix 
elements {N\qq\N), which are related to the change in the nucleon mass due to non-zero 
quark masses, denoted by cxo, and therefore to the vr-nucleon a term, S^^r as discussed later, 
and the axial-current matrix elements (A^|^7/^755'|A^), which are related to the quantities 
As^-* , a^^^ and a^g^ , as also discussed later. We find that the uncertainties in the elastic 
scattering cross section induced by the uncertainties in the quark masses, apart from the 
top quark, are negligible. However, cross section uncertainties induced by the uncertainties 
in the matrix elements {N\qq\N) and ( ^7/^75 i?|A^) are very important, as we discuss below. 
In particular, the uncertainties induced by our ignorance of S,rjv are particularly important. 

For illustration, we work within the minimal supersymmetric extension of the Standard 
Model (MSSM) with conserved R parity, and assume that the astrophysical cold dark matter 



is provided by the lightest neutralino x uM- We further assume a constrained MSSM 



(CMSSM) framework, in which the soft supersymmetry-breaking parameters mi/2,'^0 and 
Aq are assumed to be universal at the GUT input scale, and restrict our attention to scenarios 
with the Higgs mixing parameter /i > and specific values of the ratio of supersymmetric 
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TABLE I: Hadronic parameters used to determine neutralino-nucleon scattering cross-sections, 
with estimates of their experimental uncertainties. 



Higgs v.e.v.s tan/3 [14l| . We illustrate our observations by studies of some specific CMSSM 



benchmark scenarios [l5j, and also by surveys along strips in the {mi/2, mo) plane for tan f3 = 
10, 50 along which f—x coannihilation maintains the relic neutralino density within the range 
favoured by WMAP and other experiments . 



We find that the spin- independent cross section may vary by almost an order of magnitude 
for 48 MeV < St^jv < 80 MeV, the ±2-0" range according to the uncertainties in Table [H This 
uncertainty is already impacting the interpretations of experimental searches for cold dark 
matter. Propagating the ±2-0" uncertainties in A^'^'*, the next most important parameter, 
we find a variation by a factor ~ 2 in the spin-dependent cross section. Since the spin- 
independent cross section may now be on the verge of detectability in certain models, and 
the uncertainty in the cross section is far greater, we appeal for a greater, dedicated effort 
to reduce the experimental uncertainty in the n-nucleon a term Tit,n. This quantity is not 
just an object of curiosity for those interested in the structure of the nucleon and non- 
perturbative strong-interaction effects: it may also be key to understanding new physics 
beyond the Standard Model. 



II. SUPERSYMMETRIC FRAMEWORK 



We briefly review in this Section the theoretical framework we use in the context of the 
MSSM; for more comprehensive reviews, see, e.g. 



17, 18 



The neutralino LSP is the lowest- 
mass eigenstate combination of the Bino B, Wino W and Higgsinos Hi 2, whose mass matrix 
is diagonalized by a matrix Z: diag(m^^ ^) = Z*NZ^^. The composition of the lightest 
neutralino may be written as 



X = Z^iB + Z,2W + Z 



'X3 



Hi + Z^4 



H,. 



(1) 



As already mentioned, we work here in the context of the CMSSM and assume universal- 
ity at the supersymmetric GUT scale for the gaugino masses, mi/2, as well as for the soft 
scalar masses, mn, and tri-linear terms, Aq. Our treatment of the sfermion mass matri- 
ces M follows [l9|, 120]. The sfermion mass-squared matrix is diagonalized by a matrix 77: 
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diag(m1, rn^) = rjM'^r] ^, which can be parameterized for each flavour / by an angle 6f. We 
ignore here all possible CP-violating phases. The diagonalization matrix can be written as 



cos6'/ sin 6'/e*'''-' 
sin 6'/e~*'''-'' cos 6*/ 



) - 




Vl2 \ 






V22 ) 



(2) 



The magnitudes of /i and the pseudoscalar Higgs mass are calculated from the elec- 
troweak vacuum conditions using mt = 171.4 GeV, except as noted in Section HVl 

The only four-fermi Lagrangian contributions for describing elastic x-nucleon scattering 
obtained from the MSSM Lagrangian which are not velocity dependent, and hence relevant 



for relic dark matter scattering, are 191 ]: 



= a2iXYl^X(lilf^l^(lt + asiXXQiQi- (3) 

This Lagrangian is to be summed over the quark generations, and the subscript i labels 
up-type quarks {i = 1) and down- type quarks {i = 2). The coefficients are given by: 



^2 



4m?, cos^ 9w 



Z,f-\Z,f]^ (4) 



and 



+ D.Q (-^ - -i-) Re Iki (gZ^2 ' s'ZxO] I ■ (8) 

\ ''''H2 ''''Hi J ) 



where 



= r]^i— ^ r^^^agZ^ 

2mwDi ^ 



V. - Vlr + 9TuZ,2) + ^22^^^, (6) 

where Ui.T^i denote hypercharge and isospin, and 

^li = Z^z{Zxi)-, ^2i = Zy^i{—Z^^), (7) 

Bi = sin/3(cos/?), Cj = sinQ;(cosQ;), Z^j = cos a(— sin a) (8) 

for up (down) type quarks. We denote by 711^2 < the masses of the two neutral scalar 
Higgs bosons, and a denotes the neutral Higgs boson mixing angle. 
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III. HADRONIC MATRIX ELEMENTS 

The elastic cross section for neutralino scattering off a nucleus can be decomposed into a 
scalar (spin-independent) part obtained from the a^i term in Eqn. ([3]), and a spin- dependent 
part obtained from the a2i term. Each of these can be written in terms of the cross sections 
for elastic scattering for scattering off individual nucleons, as we now review and re-evaluate. 

A. Spin-Independent Term 

The scalar, or spin- independent (SI), part of the cross section can be written as ^ 

asi = ^[Z/p + (A-Z)/„f , (9) 

71 

where is the x-nuclear reduced mass and 

A = y fi''^^ + -dV y ^ (10) 

g=u,d,s ^ g=c,b,t ^ 

for = p or n. The parameters /^^^ are defined by 

m^4f = {N\m^qq\N) = m^B[''\ (11) 



where [22, [23 



/S' = 1 - E A?- (12) 

5=u,d,s 

We take the ratios of the light quark masses from 



Tfl Tfl 

— = 0.553 ±0.043, — = 18.9 ±0.8. (13) 
md md 



We take the other quark masses from |5[], except for the top mass, which is taken from the 
combined CDF and DO result [6]. These masses, as well as other experimental quantities 
that will arise in the calculation of the hadronic matrix elements, appear in Table [B 



Following [2J], we introduce the quantity: 



r(p) r(p) 
J^ip) _ ^(P) 



1.49, (14) 



which has an experimental error that is negligible compared with others discussed below, 
and the strange scalar density 



^ This expression is valid in the zero- momentum-transfer limit. For non-zero momentum exchan ge, the 



expression must include a form factor due to the finite size of the nucleus. See, for example, Ref. 21 1. 
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In terms of these, one may write 

B!t^ _ 2 + ((^-l)xy) 



# 2xz-iiz-l)xy)- ^^^^ 
Proton and neutron scalar matrix elements are related by an interchange of -Bu and B^, i.e., 

B^:'^=bI^\ 5^ = S(p), and Bi'^^ = Bi^\ (17) 
The TT-nucleon sigma term, S^jv, may be written as 

= ^(m, + ma) x (sf ) + )) , (18) 

and the coefficients /t, may be written in the forms: 

" ^ " -^(l+^)(l+f) ' ^ ^ 

" ^ " --(l + ^)(l + fe) ' ^ ^ 

~ ~ m^(l + ^) ' ^ ^ 

where we have dropped the (A^) superscript from /j-^ and Bg. 

The effect of the uncertainties in the /ju were considered in 25, 2^ and we were motivated 



to reconsider [9| the value of y by recent re-evaluations of the vr-nucleon sigma term Y^^n, 
which is related to the strange scalar density in the nucleon by 

y = l- ao/S,^ . (22) 

The value for ao given in Table [II is estimated on the basis of octet baryon mass differences 
to be (Jo = 36 ±7 MeV 0,117, 28, 29]. Recent determinations of S^jv have found the following 



values at the Cheng-Dashen point t = +2ml {8l]: 

EcD = (88 ± 15, 71 ± 9, 79 ± 7, 85 ± 5) MeV . (23) 

These should be corrected by an amount —Ar — ~ —15 MeV to obtain S,rjv- Assuming 
for definiteness the value T^cd = 79 ± 7 MeV, one finds 

S^^ = (64 ± 8) MeV . (24) 

This is the range generally considered in this paper, though one could even argue for a larger 
uncertainty, and we also discuss the implications if S^^v = <^o, i.e., y = 0. 

B. Spin-Dependent Term 

The spin-dependent (SD) part of the elastic x-nucleus cross section can be written as ^ 

asD = -GlmlAMi.J+l), (25) 

TT 



^ As with the SI cross section, this expression apphes in the zero momentum transfer hmit and requires 
an additional form factor for finite momentum transfer. This form factor may have a small but non-zero 
dependence on Cp and an- 
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where rrir is again the reduced neutrahno mass, J is the spin of the nucleus, 

A = i(ap(Sp) + a,(^,)) , (26) 



and 

The factors Ag^^ parametrize the quark spin content of the nucleon and are only significant 
for the light (u,d,s) quarks. A combination of experimental and theoretical results tightly 
constrain the linear combinations [sl 

^(P) ^ ^(p) _ ^(P) _ ^ 2695 ± 0.0029 (28) 

and 0, [n 



a 



= A(,p) + Af> - 2A(P) = 0.585 ± 0.025. (29) 

However, the individual A,^"* are relatively poorly constrained; using the recent COMPASS 
result I4I , 

A(p) = -0.09 ± 0.01 (stat.) ± 0.02 (syst.) ^ -0.09 ± 0.03 , (30) 

where we have conservatively combined the statistical and systematic uncertainties, we may 
express A^^^ as follows in terms of known quantities: 

A(f) = ^ (4^) + ) + A(P) = 0.84 ± 0.03 , (31) 
A^P) = i (4P) - a^P^) + A(P) = -0.43 ± 0.03 . (32) 

The above two uncertainties and that of A^^^ , however, are correlated and we shall instead 
use the independent quantities a^^^ and a^^ when appropriate. These values differ by ap- 
proximately 2a from those used in [2^ and we will explore the impact of this change on the 
magnitude and ratios of the spin-dependent cross sections. The proton and neutron scalar 
matrix elements are related by an interchange of Au and Ad, or 

A('^) = Af, A(")=A(p), and A^'^) = A^p). (33) 



IV. SCATTERING CROSS SECTIONS, RATIOS AND UNCERTAINTIES 



Direct and some indirect dark matter detection techniques involve the scattering of a 



WIMP off a nucleus. Direct detection experiments, such as CDMS |30|,l31i, XENONIO 



331], ZEPLIN-II |34j . |35|], and KIMS [361] aim to detect dark matter via scattering of relic 



Kamiokande [38|] and AMANDA/IceCube 



neutralinos off nuclei inside the detectors [37]. Indirect detection experiments such as Super- 



on the other hand, search for high-en ergy 

13 or Earth 0. 



neutrinos produced in WIMP annihilations at the center of the Sun 
The annihilations are the result of WIMPs accumulating at the centers of these massive 
bodies due to galactic WIMPs scattering off nuclei in these bodies and losing enough energy 
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to become gravitationally bound. Subsequent scatters then cause the WIMPs to fall to 
the cores, where the local density is enhanced. Other indirect detection methods search for 
WIMPs that annihilate in the Galactic Halo or near the Galactic Center, where they produce 



neutrinos, positrons, or antiprotons that may be seen in detectors on the Earth [42|, 143 
However, such annihilations do not involve scattering off nuclei, and will not be discussed 
here. 

The rate of such scattering events and the relative sensitivities of different experiments 
are dependent upon the four x-nucleon scattering cross-sections (spin-independent (SI) and 
spin-dependent (SD) cross-sections for each of the proton and neutron). The interpretation 
of any positive signal in such an experiment requires an understanding of the precision 
that such a signal can be correlated with underlying parameters within a given theoretical 
framework {e.g., mg, mi/2, and tan/5 in the CMSSM). 

Direct and indirect detection signals are inherently proportional to pqct^jsi, rather than 
Po or a^iy independently, where po is the local density of dark matter (assumed here to be 
dominated by relic neutralinos) and cr^jv is any of the four x-nucleon cross sections. Since 
the local dark matter density has not been directly measured, it is typically inferred from 
galactic dynamics and A^-body simulations. By convention, experimental results are often 
presented for a neutralino density of 0.3 GeV/cm^, with the implicit understanding that the 
following significant uncertainties exist in this value. In the case of a smooth distribution 



of galactic dark matter, po is estimated to be 0.2-0.4 GeV/cm^ [441 . |45| for a spherical halo, 
but it may be somewhat higher, up to 0.7 GeV/cm^ 46|, |47|], for an elliptical halo; see 



Ref. [45| for a discussion of the difficulties in determining this value. Models of the galaxy 
based upon hierarchical formation [is!], which do not assume a strictly smooth distribution 
of dark matter as do the above estimates, suggest that the local density may be as low as 
0.04 GeV/cm'^, but that 0.2 GeV/cm^ is a more reasonable lower limit. With this hierachical 



formation, the halo may contain substructure such as clumps [49|, |50| or tidal streams from 
galaxy mergers 51| that can increase the dark matter density. An analysis of the Sagittarius 
dwarf galaxy, whose leading tidal tail passes through the galactic disk somewhere near our 
location, estimates it may contribute an additional 0.001-0.07 GeV/cm^ to the local dark 
matter density |52]. Without detailed knowledge of the size and prevalence of all such 
substructure (particularly, the locally present substructure), it is difficult to place upper 
limits on the local dark matter density; however, an increase by greater than a factor of a 
few over the above estimates are unlikely j48|. On the other hand, one should in principle 
allow for the possibility that there may be some additional source of cold dark matter, such 
as axions, in which case the neutralino density would be less than po, though we do not 
consider this possibility here. Without a precise determination of po, it is difficult to place 
precise limits on the x-nucleon cross sections from direct detection experiments. However, 
the ratios of cross-sections may be unambiguously determined, so we focus on these ratios 
in the following. 

We take the fiducial values of the experimental quantities to be those obtained using 
the central values of the parameters in Table [B For illustration, cross-sections and their 
ratios are given in Table [TTl for a few benchmark models which are essentially the well- 



studied benchmark models C, L, and M 15 . All of the benchmark models lead to relic 



densities within the WMAP preferred range of flh'^ = 0.088 — 0.120 53|], and satisfy most 



phenomenological constraints ^. Points C and L are points along the coannihilation strip 



^ The exception being a possible failure to account for the discrepancy between the theoretical and experi- 
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(at low tan/3 = 10 and high tan (3 = 50, respectively) where the masses of the neutralino 
and stau are nearly degenerate. Point M is in the funnel region where the mass of the 
neutralino is roughly half that of the Higgs pseudoscalar. Figure [1] shows the cross sections 
along the WMAP-allowed strips of the {mi/2,mQ) planes for tan/3 = 10 and 50 [l|, |l6|, 56 



For tan (3 = 10, this is the coannihilation strip and includes point C, whereas for tan/3 = 50, 
the strip is formed by coannihilations at low mi/2, turning into the rapid annihilation funnel 
at larger mi/2. This strip includes points L and M. 

Current experimental upper limits on cross sections are also given in Figure [TJ All 
these limits assume a local neutralino density of 0.3 GeV/cm^, and should be rescaled by 
(0.3 GeV/cm^/po for other values of po- We display the SI limits given by CDMS and 
XENONIO |32| under the reasonable CMSSM assumption that o"^p,si ~ crxn,si; relaxing that 
assumption would weaken the given limits by at most a factor of ~5. The same assumption 
cannot be made for the SD cross-sections, and experiments typically constrain only one of 
<^xp.SD or cr^n SD, not both. The direct detection experiment KIMS [36| provides a limit for 



Cxp.SD' while CDMS [31| and XENONIO [33| (taken from [57[) constrain cr^n.sD- Since the 



Sun is primarily composed of hydrogen, the capture rate of neutralinos is particularly sensi- 
tive to ct-^p.sd; thus, indirect detection experiments are able to constrain this cross-section. 
The large mass of the Sun and the large size of the detector allows Super-Kamiokande 38 
to provide a significantly better limit on cr^p,sD than the direct detection experiments. It 
is apparent from Figure [1] that experiments are beginning to probe the CMSSM parameter 
space through SI scattering, particularly at large tan /?, but are not yet sensitive enough to 
detect SD scattering. 

From our analysis of the benchmark models and along the WMAP strip, several general 
features of the x-^ucleon cross sections are apparent for any given CMSSM model. The 
SI cross sections for x-proton (o";^p,si) and x-neutron (a^^n.si) scattering are typically very 
close, within ~5% of each other; close enough that they are virtually indistinguishable in 
Figure [H Likewise, the spin- dependent x-proton (c'"^p,sd) and x-neutron (cr^^n.so) scattering 
cross sections are similar, differing by at most a factor of 2 or 3. In general, the SD x- 
nucleon cross section is much larger than the SI one, by O{10^ — 10^) or more. However, 
we recall that the SI cross section for scattering off a nucleus, Eqn. contains a factor 
of the number of nucleons squared, whereas the SD cross section, Eqn. fl2^ . is proportional 
to the square of the spin, which does not grow with increasing nuclear mass. Consequently, 
heavy elements such as Ge and Xe are actually more sensitive to SI couplings than to SD 
couplings. 



A. Confidence Intervals 

To determine the level of uncertainty in the cross sections and ratios obtained for a given 
model, we generate confidence intervals about the central values of the resulting cross section 
and ratio distributions by Monte Carlo sampling over the parameters in Table [H assuming 
independent Gaussian errors ^. Results are given in Table HTl for benchmark models C, L, 



mental results for anomalous magnetic moment of the muon [541 . l5a | . Points C and L are consistent with 
{g — 2)^, whilst the contribution from point M is too small. 
^ The distributions for the cross sections generated by the Monte Carlo are not even approximately Gaussian 
or symmetric as is evidenced by the curvature seen in Figure^ Hence, confidence intervals are used instead 
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FIG. 1: The neutralino-nucleon scattering cross sections along the WMAP-allowed coannihila- 
tion strip for tan/3 = 10 (left panels) and coannihilation/funnel strip for tan/3 = 50 (right panels) 
using the parameters in Table [H CMSSM predictions for the spin-independent (SI) cross sections 
are shown in the upper panels (^xp.si in solid/blue, o";^n,si in solid/red), along with experimen- 
tal contraints from CDMS (dotted/black) and XENONIO (dashed/black) , where the constraints 
apply to either SI cross section under the reasonable MSSM assumption c^p.si ^ '^xn.si- Note 
that c^p.si and c^n.si are nearly indistinguishable at the scale used for these plots. The lower 
panels show CMSSM predictions for the spin-dependent (SD) cross-sections ((7^p,sD in solid/green, 
c"xn,SD in solid/magenta). Experimental constraints for o"xp,sD are shown for KIMS (dotted/green) 
and Super-Kamiokande (dashed/green), while constraints for a^n.SD are given by CDMS (dot- 
ted/magenta) and XENONIO (dashed/magenta). Limits are based upon a local neutralino density 
of 0.3 GeV/cm3. 

and M at confidence levels (CL.) of 68.3% and 95.4%. These confidence intervals take into 
account only hadronic uncertainties, not any uncertainties arising within the supersymmetric 
models themselves. We recall, in particular, that the estimates for the a2q and asq coefficients 
of the effective Lagrangian in Eqn. ([3]) contain uncertainties that may only be reduced by 
calculating radiative corrections to the effective scattering operators, work that is beyond 



of the averages and standard deviations of the generated distributions. 
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0.78 - 2.14 




0.75 - 1.80 




0.71 - 1.46 




ctxp.sd/ctxp.si 


770 




77 




187 




68.3% C.L. 


480 - 1350 




46 - 151 




121 - 319 




95.4% C.L. 


320 - 2730 




30 - 373 




83 - 616 




0'xn,SD/crxn,SI 


960 




86 




186 




68.3% C.L. 


610 - 1660 




51 - 163 




120 - 313 




95.4% C.L. 


410 - 3370 




33 - 392 




83 - 601 





TABLE IL Neutralino-nucleon scattering cross sections and ratios for benchmark models C, L, 
and M. Confidence intervals are given for confidence levels (C.L.) of 68.3% and 95.4%, using the 
hadronic parameter uncertainties in Table HI 

the scope of this paper. These confidence intervals take into account only the explicit 
dependence of nit in Eqns. (HI), (E]), & ( fTOj) . but not in MSSM terms such as 9t, rriA, and 
|yu|, which are calculated using the central value of 171.4 GeV. The confidence intervals are 
not significantly affected by this incomplete treatment of the top quark mass uncertainty, 
apart from a slight underestimate of the ct^tv.sd interval widths; we shall address below the 
full MSSM dependence on the top quark mass. 

As indicated by the benchmark models, the hadronic uncertainties induce significant 
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f^xp,si (pb) 


-^xp.SD (pb) 


0-xn,Sl/0-XP,SI 


'7xn,SD/'7xp,SD 


fiducial value 


2.85 X 10 ^ 


2.19 X 10 ° 


1.029 


1.28 




±3.5% 


~ 0." 


±0.08% 


~ 0. 




~ 0. 


~ 0. 


~ 0. 


~ 0. 


nn 1 nn 1 
'"-s/ '"'d 


+ 5 2% 

-i—O. ^ /o 


~ 


±0 07% 


~ 


rric 


~ 0. 


~ 0. 


~ 0. 


~ 0. 


mb 


±0.1% 


±0.04% 


~ 0. 


±0.01% 


mt 


±9% 


±9% 


±0.005% 


±2% 




+34% 
-27% 




+1.2% 
-0.7% 






+45% 
-32% 




+0.7% 
-0.4% 








±0.5% 




±0.06% 






±2.2% 




±4.2% 






+14% 
-12% 




+30% 
-21% 



Entries listed as ~ 0. have a relative variation less than 10 ^ that of the fiducial value. 



TABLE III: Relative uncertainties in (Txp.sIi fxp,SD, crxn,si/cxp,Si) and (7;^n,SD/o'xp,SD for Model 
C due to each of the parameters in Table HI The quoted uncertainties correspond to 68.3% C.L. 
confidence intervals relative to the fiducial values. Uncertainties in o";^n,si and o";^n,SD (not shown) 
are comparable to those in c^^-p^si and cr^p,sD5 respectively. Due to the non-linear dependence of 
the cross sections on some of the parameters over the ranges of interest, some confidence intervals 
are not symmetric about the fiducial values. 

uncertainties in the cross sections. At the 68.3% C.L., the SI cross sections vary by ~3; at 
the 95.4% C.L., they vary by an order of magnitude. Such a variation is larger than might 
naively be expected given the much smaller relative uncertainties of the parameters listed 
in Table [U the largest of which is at 5±2 MeV ^ . Clearly, such a large uncertainty in the 
cross sections would make it difficult to narrow down CMSSM parameters on the basis of 
any single experimental signal. 

Whilst the two SI cross sections can vary greatly, their ratio does not: in our Monte 
Carlo analysis, cr^n,si/c"xp,si varies by only a few percent. It is worth noting, however, that 
in some models {e.g. model L), a^^n.si may exceed cr^p,si by more than 10% at the 95.4% 
C.L., a difference that may be of experimental interest. In addition, even at the 95.4% C.L., 
c"xn,si is always larger than cr^p,si- 

The variations in the two SD cross sections are not so large: ~20% and ~40% at the 
68.3% and 95.4% C.L.'s, respectively. The ratio of the two cross sections, however, has a 
larger variation, typically by a factor of 2 or 3, due to an anti-correlation in the variations 
of the two cross sections that we discuss below. 

To demonstrate the contributions of individual parameters to the cross-section uncertain- 
ties, we have determined 68.3% C.L. confidence intervals for variations in the parameters 
treated singly, with the remaining parameters fixed at their central values, as presented in 



^ We show below that, whilst having the largest relative uncertainty, md makes a negligible contribution to 
the uncertainties in the cross sections. 
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Table UTT] for benchmark model C. The intervals are given as variations relative to the fidu- 
cial values. For several parameters, notably ctq, E^jv, and A^'''*, the variations in the cross 
sections are non-linear with respect to the parameters over the values of interest, resulting 
in confidence intervals that are not symmetric about the fiducial values. The variations due 
to the top mass here include the dependence of MSSM parameters such as 6't, Aq, and 
on rrit, which was neglected in the full confidence intervals of Table HIl 

It is clear that from Table IIIII that the uncertainties in the quark masses and their 
ratios, apart from the top mass, make almost negligible contributions to the cross-section 
uncertainties. The SI cross-section uncertainties are dominated by contributions from (Tq 
and S^jv; and the SD cross-section uncertainties mainly arise from uncertainties in As''''. We 
examine these in the following sections. 



B. Quark Masses and Mass Ratios 

The most significant quark mass dependence of the x-nucleon scattering cross sections 
is that on the top-quark mass, which induces uncertainties of ~10% in all four <j^n- These 
uncertainties arise not from the appearance of rrit in Eqns. dlj), dSD, & f|TU]) (hereafter referred 
to as the explicit rrit dependence), but in the calculations of CMSSM parameters, e.g., 
and However, variations in rut rescale all four cross- sect ions in the same manner, so 
that ratios of the cross-sections are essentially independent of rrit. The precision of the top 
mass measurement will continue to improve once the LHC begins taking data, with a 1 GeV 



uncertainty a possibility in as little as one year of low-luminosity running [58j. Thus, we 
can expect the mt-induced uncertainties in a^^ to fall from ~10% to ~5% or lower within 
a few years, well below uncertainties induced by other parameters. 

The SI cross sections are sensitive to the ratios of the light quark masses mu/mj and 
mg/md, more than to the overall scale of the light quark masses, which is fixed by the 
relatively poorly constrained value of rrid. In the CMSSM, 6f ^ for the light quarks 
{u,d,s) and the charm quark (c), so that r] in Eqn. ([2]) is nearly diagonal. For diagonal r], 
a^q is proportional to rriq, and hence asq/nig in Eqn. ffTOl) is independent of the quark mass. 
Any dependence on the light quark masses must come through the f^^^^ terms. However, 
by Eqns. f lTI?]) - fl2T]) . those terms depend only on mass ratios. Since there is no dependence 
on md (for fixed mass ratios), its relatively large uncertainty does not translate into any 
significant uncertainties in the SI cross sections. 

Whilst the mass ratios induce uncertainties of a few % in the cross-sections, the induced 
uncertainties in the ratio cTj^n,si/o"xp.si are only ~0.1%. This is a consequence of the fact that 
the predominant contribution to Eqn. (fTOl) comes from the strange-quark term f^^^^cyss/'ms, 
which is common to both /p and f^. Thus, variations in the the quark mass ratios induce, 
via Eqn. (pTi) . comparable (and correlated) variations in both cT^p,si and cr^n.si, resulting in 
only small variations in the ratio of those two cross sections. 

The r] matrices for the bottom and top quarks, unlike the lighter quarks, are not ap- 
proximately diagonal, so that asq/rrig does have a quark mass dependence. However, the 
precisions to which these quark masses are known lead to only small uncertainties (<0.1%) 
in the SI cross-sections (here, we refer only to the explicit rrit dependence in these terms, 
not to the MSSM parameter dependence discussed previously). 

The only quark mass dependence of the SD cross sections in Eqn. ( 1271) comes from the 
light quark masses in a2q {q = u,d,s). However, mass-dependent terms are suppressed by 
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Model 


C 


L 


M 


= 36 MeV: 








o-xp.si (Pb) 


3.40 xlO-i° 


1.38 xlO-9 


1.78 xlO-ii 


o"xn,si (pb) 


3.67 xlO-i° 


1.61 xlO~^ 


1.89 xlO"^^ 


0-xn,Sl/o-xP,SI 


1.080 


1.170 


1.065 


= 45 MeV: 








c^xP.si (pb) 


8.80 xlO-i° 


5.55 xlO-9 


4.23 X 10-11 


CTxn.Si (pb) 


9.24 xlO-i° 


6.02 xlO-9 


4.41 X 10-11 


0-xn,Sl/0-xP,SI 


1.050 


1.085 


1.043 


E^^r = 56 MeV: 








o-xp.si (pb) 


1.88 xlO-9 


1.45 xlO-s 


8.64 X 10-11 


o-xn,si (pb) 


1.95 xlO-9 


1.52 xlO-s 


8.91 X 10-11 


0-xn,Sl/0-xP,SI 


1.035 


1.053 


1.031 


= 64 MeV: 








o-xp,si (pb) 


2.85 xlO-9 


2.36 xlO"*^ 


1.28 xlO-i° 


o-xn,si (pb) 


2.93 xlO"^ 


2.46 xlO"*^ 


1.32 xlO-i° 


0-Xn,Sl/o-xp,SI 


1.029 


1.042 


1.026 


S^iv = 72 MeV: 








c^xP,si (pb) 


4.01 xlO"^ 


3.49 xlO"*^ 


1.78 xlO-i° 


0-xn,SI (pb) 


4.11 xlO"^ 


3.61 xlO-® 


1.82 xlO-i° 


C^Xn.Sl/o-xP-SI 


1.025 


1.035 


1.022 


= 84 MeV: 








o-xp.si (pb) 


6.13 xlO"^ 


5.61 xlO"*^ 


2.69 xlO-i° 


o-xn.si (pb) 


6.26 xlO"^ 


5.76 xlO"*^ 


2.74 xlO-i° 


0'xn,Sl/0"xp,SI 


1.021 


1.028 


1.019 



TABLE IV: Spin-independent neutralino-nucleon scattering cross sections in the benchmark mod- 
els for several values of S^jv- 

a factor of 0{mq/mw) {0{mg/m'^) for diagonal r] matrices) relative to mass-independent 
terms arising from the Xi, Yi, Wi, and Vi factors given by Eqn. ([6]). Thus, the SD cross 
sections are nearly independent of the quark masses. 

C. Spin-Independent Parameters and Cross Sections 

The greatest impediment to an accurate determination of the SI cross sections for any 
given MSSM model comes from the Sttjv and cxo parameters. As shown for benchmark model 
C in Table llllt each of these two parameters induces uncertainties of ~30% or more in cr^p^si 
at the 68.3% C.L.; cr^n.si (not shown) has similar induced uncertainties. The large confidence 
intervals for the SI cross sections in Table [III are almost entirely due to the uncertainties in 
these two parameters. 

We focus the discussion here mainly on Sttjv rather than ctq as the S-n-jv — 64 ih 8 MeV 
result is significantly larger than previous estimates for that parameter. Indeed the range of 
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rN 



(MeV) 



FIG. 2: The spin- independent neutralino-nucleon scattering cross section as a function of ^-j^n 
for benchmark models C, L, and M. Note that Cxp.si and o";^n,si are nearly indistinguishable at the 
scale used in this plot. 
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FIG. 3: The spin-independent neutralino-nucleon scattering cross section ratio along the WMAP- 
allowed strips for tan/3 = 10 (left panel) and tan/3 = 50 (right panel) for several values of S^jv- 
Experimental constraints from CDMS and XENONIO are also shown. 



estimates of the central value for S^^jv is far greater than the typically quoted uncertainty. 
In view of this, we also include below some results for lower values of S^jv- 

In Figure O we show the S^jv dependence of (J^n,si for the benchmark models, and 
Table HVl gives the cr^N,si values for those models for selected values of T,.,,j^. All the other 
parameters are set at their fiducial values (Table [11). From the minimal value for S^jv (ctq = 
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36 MeV) to the 2-a upper bound (80 MeV), ct^tv.si varies by more than a factor of 10 (as 
much as a factor of 35 for model L). At these benchmarks and in other models of interest, for 
larger values of S^riv (X'-kn 7^ ctq), the majority of the contribution to /p in Eqn. (fTOj) comes 
from the strange quark term, with f!j^^ (x y oc S^jv — o"o; so that o"^p,si ~ C^-kn — <7o)^- Thus, 
the SI cross sections are particularly sensitive not just to and Uo, but to their difference. 
For smaller values of Et^jv {^nN ~ o'q), the strange contribution no longer dominates, but a 
strong dependence of cr^Ar,si on S^jv does remain. 

In Figure [31 the SI cross sections are shown along the WMAP allowed coannihilation strip 
for tan/5 = 10 and coannihilation/funnel strip for tan/3 = 50 for the Sttjv reference values of 
36 MeV (no strange scalar contribution), 64 MeV (central value), and 80 MeV {2-a upper 
bound); CDMS and XENONIO limits are also given. As with the benchmark models, a 
factor of ~10 variation occurs in cr-^Ar,si over these Sj^iv reference values for any given model 
along the WMAP strip. 

Such large variations present difficulties in using any upper limit or possible future pre- 
cision measurement of cr^N,si from a direct detection signal to constrain the CMSSM pa- 
rameters. The present CDMS and XENONIO upper limits have (almost) no impact on the 
WMAP strip for tan/? = 10(50), if one makes the very conservative assumption that ctq = 
36 MeV {y = 0). On the other hand, mi/2 ~ 200 GeV would be excluded for tan/? = 10 if 
S^^ = 64 or 80 MeV. This excluded region would extend to mi/2 ~ 300 GeV for tan/? = 50 
if S^riv = 64 or 80 MeV. Thus, the experimental uncertainty in E^jv is already impinging on 
the ability of the present CDMS and XENONIO results to constrain the CMSSM parameter 
space. 

Looking to the future, a conjectural future measurement of cr^^p.si = 4 x 10^^ pb would 
only constrain mi/2 to the range 600 GeV < mi/2 < 925 GeV if tan/5 = 10 and 1100 GeV 
< mi/2 < 1400 GeV if tan/? = 50, for the 1 -0" Sjrjv range of 64 i 8 MeV. If smaller values 
of Ejrjv are also considered, down to ctq — 36 MeV [y — 0), these constraints would weaken 
to 350 GeV < mi/2 < 925 GeV and 550 GeV < mi/2 < 1400 GeV for tan ;3 = 10 and 50, 
respectively, ruling out only the smallest values of mi/2 ^■ 

As the SD cross sections are independent of S^^, the SD/SI cross-section ratio exhibits 
the same behavior with respect to S^^ as the SI cross-section does (albeit inverted), as shown 
in Figure HI The ratio is also shown along the corresponding WMAP strips in Figure [5] for 
several values of Sttat ^- For small S,rjv, the ratio is large, almost 10^ for model C and 2000 
for models L and M at Sttat = 36 MeV, suggesting that spin-sensitive experiments have a 
better chance of detecting neutralinos than spin-insensitive experiments. At the 2-a upper 
limit on S^^, however, the ratio falls to 500 for model C, 50 for model L, and 120 for model 
M. Recalling that the SI x-nuclear scattering cross section scales as the number of nucleons 
squared, SI scattering may actually dominate over SD scattering for heavier nuclei in the 
cases of these lower SD/SI cross-section ratios. 

The dependence of the SI x-^/x-P cross-section ratio on is shown in Figure El for 
the benchmark models, with ratios at selected S^^r given in Table [IVl The cr^n,si/cxp,si 
ratio is shown along the WMAP coannihilation strips in Figure [7] for several values of T^t^,^. 



^ Moreover, as noted previously, however, detection signals only measure poa^N, so a^N can only be deter- 
mined from a signal to the precision that the local dark matter density is known. 

^ The non-monotonic dependence on mi/2 seen in the right panel of Figure [5] for tan/? — 50 is due to 
an enhancement of the Higgsino components, Z-^^ and Z^^, at low mi/2- As a result, the Z-exchange 
contribution to the spin-dependent cross section increases rapidly. 
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FIG. 4: The ratios of the spin-dependent and spin-independent neutrahno-nucleon scattering 
cross sections as functions of S,rjv for the benchmark models C, L, and M. 
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FIG. 5: The ratios of the spin-dependent and spin- independent neutrahno-nucleon scattering 
cross sections along the WMAP-allowed strips for tan (3 = 10 (left panel) and tan /3 = 50 (right 
panel) for several values of S^jv- 



Whilst cT^p,si and cr^^n.si each depend significantly on the value of Sttiv, their ratio has only 
a mild E^jv dependence. As noted earlier in this section, the strange-quark term dominates 
in Eqn. flTU]) for S^riv 7^ ctq and, since this term contributes identically to both /p and /n 
(neglecting nucleon mass differences), cr^n,si/cxp,si ~ 1; independent of E^jv- For E^^v ~ o"o, 
where the strange contribution is no longer significant, we generally find that Jt^, fx^ <^ 1 
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FIG. 6: The ratios of the spin-independent neutralino-neutron and neutralino-proton scattering 
cross sections as functions of S^jv for benchmark models C, L, and M. 
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FIG. 7: The ratios of the spin- independent neutralino-neutron and neutralino-proton scattering 
cross sections along the WMAP-allowed strips for tan (3 = 10 (left panel) and tan /3 = 50 (right 
panel) for several values of S^jv- 



so that fxG ~ 1 via Eqn. f|T2l) and the right (heavy-quark) summation term in Eqn. ffTOj) 
dominates over the left (light-quark) summation term. Since the heavy-quark terms are 
identical in /p and /n, we again have o"^n,si/c"xp,si ~ ^■ 

The f!j^^ and f^'^ terms are small but not entirely negligible, and the fact that f!jf^ 7^ 
for g=u,d implies that /p 7^ /n in general. The result is a small shift in the SI x-ii/X"P 
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cross-section ratio away from unity. The shift becomes greater as S^jv ~^ '^o, since the 
dominant strange-quark contribution disappears and allows the proton/neutron asymmetry 
in the smaller up and down quark terms to become relevant. For the benchmark models in 
Figure E] and along the WMAP coannihilation strips in Figure [TJ cr^n.si is ~5% larger than 
o";,.p,si in the 2-a St^jv range 64 ± 16 MeV. For E^jv ~ ctq, however, the difference can become 
as large as 15 to 20 %. 

It can be shown that o";^n,si/c"xp.si > ^ ^'^^ \(^3d\ > l^sul- Since oc = 1/ sin/? 

and asd oc 1/-Bd = 1/ cos/3, that condition is generally satisfied for the tan (3 ^ 1 models 
examined here. Hence, the neutron always has a slightly larger SI cross section than the 
proton in these models. 

When including all the uncertainties, as shown by the Monte Carlo results in Table [Tll 
a small but non-zero asymmetry in the proton and neutron SI cross-sections is expected 
to occur at the ~ few % level. Such an asymmetry may be observable, but a problem 
arises in the current generation of experiments. To extract the relative strength of the /p 
and /n couplings in Eqn. a signal must be seen in two detectors with different ratios 
of Z and A — Z. However, to boost the scattering cross-section and, therefore, neutralino 
detection likelihood, many current experiments use heavy elements such as Ge and Xe that 
have similar ratios A — Z ^ lAZ. It will be necessary to see a signal in experiments using 
lighter elements nearer to A — Z ^ Z {i.e. Z <17) in order to determine cr;^n,si/<7xp,si- Such 
experiments are possible, but their SI sensitivities typically lag far behind those with heavier 
elements. 



D. Spin-Dependent Parameters and Cross Sections 

The determination of the SD cross sections for a given MSSM model depends on the three 
parameters specifying the spin content in a nucleon: a^^\ a^g\ and Ai^\ As demonstrated 
for benchmark model C in Table llllt uncertainties in 03'''' and a^^^ induce only < 1% and 
~ few % uncertainties in <j^n,sd, respectively, at the 1-a level. Uncertainties in the strange 
spin contribution As^\ on the other hand, induce 10 to 15 % uncertainties in (T^at^sd, and 
uncertainties in this parameter account for essentially all the width of the SD confidence 
intervals in Table [Tll Since uncertainties induced by 03'''' and a^g^ are negligible, we ignore 
these terms and focus on Ai^\ 

We give in Table |V] the SD cross sections and their ratios for the benchmark models for 
several A^^^ values ranging from -0.15 (the 2-cr lower bound) to 0.0 (no strange contribution 
to the nucleon spin). Over the 2-cr range —0.09 ±0.06, both cr^p,sD and cr^n,SD vary by ~40% 
in each of the three models, a significant variation but not so large as that induced in the 
SI cross sections by T,T,f^. However, unlike S^jv and the SI cross sections, the SD proton and 
neutron cross-sections are antz-correlated with A^'''': as the value of A^'^'' increases, o"^p,sD 
decreases whilst o"^n,SD increases. 

Figure [8] and Figure [9] show the SD cross sections and the SD/SI cross-section ratios along 
the WMAP coannihilation strips for several values of A^'''*. These figures again demonstrate 
the anti-correlation of the neutron/proton cross sections with As'''' as (for decreasing Ag^^) 

o"xp.SD and cr^^n.sD approach each other and, for low enough As^"*, cross over (o"^p,sD becomes 
the larger SD cross-section). 

Because of the anti-correlated behaviors of the SD cross sections, the SD x-n/x-p cross- 
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Model 



-0.15: 



o-xp,SD (pb) 
o-xn,SD (pb) 

C^Xn.SD/cTxp.SD 

A^P^ = -0.12: 
o-xp,SD (pb) 

f^Xn.SD (pb) 
0'xn,SD/o"xp,SD 

A^P^ = -0.09: 
o"xp,SD (pb) 

0-xn,SD (pb) 
fxn,SD/o"xp,SD 

A^P^ = -0.06: 
o"xp,SD (pb) 

0-xn,SD (pb) 

a 

A^P) = 



Xn.SD/^xp.SD 



-0.03: 

o-xp.SD (pb) 

0"xn,SD (pb) 

a 

a(p) = 



Xn,SD/crxp,SD 



0.00: 

o"xp,SD (pb) 

0-xn,SD (pb) 
Cxn,SD/o"xp,SD 



C 

2.80 xlO"*^ 
2.19 xlO-6 

0.78 

2.48 xl0"6 

2.49 xlO-6 
1.00 

2.19 xlO-6 

2.81 xlO^s 
1.28 

1.91 xlO-6 
3.14 xlO-6 
1.64 

1.65 xlO-6 
3.49 xlO-6 
2.11 

1.41 xlO-6 
3.86 xlO-6 
2.74 



2.25 xlO-6 
1.69 xlO-6 

0.75 

2.03 xl0~6 
1.89 xlO-6 
0.93 

1.82 xlO-6 
2.10 xl0"6 
1.15 

1.63 xlO-6 
2.33 xl0~6 
1.43 

1.44 xlO-6 
2.57 xlO-6 
1.78 

1.26 xlO-6 
2.81 xl0"6 

2.23 



M 

2.86 xlO"^ 
2.03 xlO-s 
0.71 

2.63 xlO"^ 
2.23 xlO-8 
0.85 

2.40 xlO-s 
2.45 xlO-s 
1.02 

2.19 xlO-s 
2.67 xlO-s 
1.22 

1.99 xlO"^ 
2.90 xlO^s 
1.46 

1.80 xlO-s 
3.15 xlO"^ 
1.75 



TABLE V: The spin-dependent neutralino-nucleon scattering cross sections in the benchmark 
models for several values of is!^\ 

section ratio is particularly sensitive to the value of As . This ratio is shown along the 
WMAP coannihilation strips in Figure [T0|, where it is apparent that the ratio varies by a 
factor of 2 to 3 at the ~ 2-a level for Ai^"*. That factor of 2 to 3 in the ratio can also be 
seen in the confidence intervals at the 95.4 % C.L., given in Table HTl 

The large uncertainty in cr^n.sD/o'xp.sD induced by /slf^ for a given model is unfortunate 
since this ratio may be one of the easiest to determine experimentally. This is because many 
elements have a spin-odd proton group ((5'p) 0) and a spin-even neutron group ((5'p) ~ 0) 
or vice versa, in which case A oc Cp or A oc On in Eqn. (l25l) . By using an odd-even element in 
a detector, any signal could essentially be entirely attributed to Op and therefore yield o"xp,sD- 
Likewise, when using an even-odd element, any signal could essentially be entirely attributed 
to On and therefore yield cr^n.sD- Even though the measurements would probably actually 
be of PoO'xp.SD and PoC"xn,SD and not of cr^p,sD and cT^n,SD alone, the ratio cr^n,SD/o"xp,SD is still 
an unambiguous and straightforward experimental measurement. The large uncertainties 
induced by A^^^ however, make it difficult to use such a signal to constrain the CMSSM 
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FIG. 8: The ratios of the spin-dependent neutrahno-nucleon scattering cross sections along the 
WMAP-allowed strips for tan (3 = 10 (left panel) and tan /? = 50 (right panel) for several values of 




FIG. 9: The ratios of the spin-dependent and spin- independent neutralino-nucleon scattering 
cross sections along the WMAP-allowed strips for tan/3 = 10 (left panel) and tan/3 = 50 (right 
panel) for several values of 



parameter space. 



V. SUMMARY AND DISCUSSION 



We have analyzed in this paper the principal hadronic uncertainties in the spin- 
independent (SI) and spin-dependent (SD) cross sections for supersymmetric relic scattering 
on protons and neutrons, using three benchmark points and two coannihilation strips as il- 
lustrations. We have found that the principal hadronic uncertainty in the SD cross sections 
is due to our lack of knowledge of the 7r-nucleon a term. In comparison, uncertainties in 
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FIG. 10: The ratios of the spin-dependent neutralino-neutron and neutralino-proton scattering 
cross sections along the WMAP-allowed strips for tan/3 = 10 (left panel) and tan/3 = 50 (right 
panel) for several values of 



quark masses and their ratios are much less important. In the case of the SD cross sections, 
the dominant uncertainty is due to our ignorance of the strange-quark contribution to the 
nucleon spin, though this uncertainty is relatively less important than that induced in the 
SI cross sections by the 7r-nucleon a term. 

This uncertainty in the vr-nucleon a term clouds very significantly the interpretation of 
searches for (and eventually measurements of) dark matter scattering on nuclei, preventing 
precise answers to the key questions: How do present unsuccessful searches constrain the su- 
persymmetric model parameter space? How accurately could a possible future measurement 
be used to refine the model parameters? This hadronic uncertainty is much larger than that 
generates by uncertainties in supersymmetric model calculations of the effective LSP-quark 
interactions, and also much larger than the astrophysical uncertainty in the local cold dark 
matter density. 

One of the great hopes in supersymmetric phenomenology is that one will eventually be 
able to use measurements at accelerators such as the LHC and/or a linear e~^e~ collider to 
calculate the relic LSP density in the Universe, and the rates for dark matter scattering. The 
hadronic uncertainty in the latter that is induced by our ignorance of the vr-nucleon a term 
limits severely the prospects for completing the second part of this programme. Specifically, 
this uncertainty is much larger than the uncertainty in calculating the relic LSP density that 
could be expected from LHC measurements in at least one benchmark model. 

We therefore plead for an experimental campaign to determine better the TT-nucleon a 
term. This quantity is certainly interesting and important in its own right and as a measure 
of the importance of strange quarks in the nucleon. However, as argued in this paper, it 
is potentially also a key ingredient in the effort to understand one of the most important 
aspects of possible new physics beyond the Standard Model. 
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